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ABSTRACT 

We report a search for z=7 Lya emitters (LAEs) using a custom-made 
Narrow-Band filter, centered at 9755A with the instrument Suprime-Cam in- 
stalled at the Subaru telescope. We observed two different fields and obtained 
two sample of 7 Lya emitters of which 4 are robust in each field. We are 
covering the luminosity range of 9.10 42 — 2.10 43 ergs _1 in comoving volumes of 
~ 4 x 10 5 and 4.3 x W 5 Mpc 3 . 

From this result, we derived possible z~7 Lya luminosity functions for the full 
samples and for a subsample of 4 objects in each field. We do not observe, in each 
case, any strong evolution between the z=6.5 and z~7 Lya luminosity functions. 
Spectroscopic confirmation for these candidate samples is required to establish a 
definitive measure of the luminosity function at z~7. 

Subject headings: cosmology: early universe, galaxies: luminosity function, mass 
function, galaxies: distances and redshifts 



1. Introduction 

Galaxies formed at high redshifts play a key role in understanding how and when the 
reionization of the universe took place. They also help constrain the physical mechanisms 
that drove the formation of the first stars and galaxies in the universe. Over the last decade, 
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the limits of the observable universe have been pushed at z ~ 6, which corresponds to ~ 90% 
the age of the universe. Their detection will allow us to probe the era of reionization. 
Selecting galaxies with strong emission lines, such as Lya line, allows us to probe the high- 
redshift Lya luminosity function. 

The search for the redshifted Lya emission at the longest possible wavelength is complicated 
by the presence of OH emission lines within the terrestrial atmosphere. This strong emis- 
sion line is responsible for the faintness limit at which celestial objects can be detected with 
ground-based telescopes at near-infrared (IR) wavelengths. Fortunately, there are spectral 
intervals with lower OH-background that allow for a fainter detection limit from the ground- 
based observations. This is known as the narrow-band (NB) imaging technique. 
This is one of the most successful methods to detect strong Lya emission lines of galaxies, 
since it relies on a specific redshift interval as well as a selected low-sky background window. 
This filter allows us for maximum detection of light from the Lya emitters at the central 
wavelength, while minimizing the adverse influences of sky emission. 

Over hundred of z > 6 LAEs have been p h otometrically sele c ted an d spectroscopically iden- 
tified in this way since iHu et al.l ( 120021 ). iKashikawa et al.l (120111 ) obtained an exhaustive 
sample of LAEs, (45 at z=6.5 and 54 at z=5.7 spectroscopically confirmed), from which a 
first robust estimate of the Lya LF was derived. This LF shows an apparent deficit compared 
to the z=5.7 Lya LF of IShimasaku et al.l ( 120061 ). corresponding to a possible luminosity evo- 
lution from z=5.7 to z=6.5 of L* z=6 6 ~ 0.4 — 0.6L* =5 7 . They conclude that the reionization 
of the unive r se is n ot been completed at 6.5. 



Ouchi et al.l (120101 ) have now obtained the largest sample to date of 207 LAEs at z=6.6 with 



the NB imaging technique. Their derived z=6.6 apparent Lya LF indicate a decrease from 
5.7 at 90% confidence level with a more dom inant decrease of l umino sity evolution (L*) than 
number evolution ($*), in agreement with IKashikawa et al.l ( 120061 ). They claim therefore 
that the hydrog en in the IGM is not highly neutral at z=6.6. 



Hu et al.l (120 lOl ) have obtained 88 z~5.7 and 30 z~6.5 Lya emitters. Their results on the 



evolut ion o f the Lya LF ar e in ag reement with previous works from iMalhotra fc Rhoads 
d2004h and IKashikawa etaD feood ). 

Ive et dJ ( B006h have first confirmed spectroscopically a z=6.96 LAE. From this result, 



Ota et al.l (120081 ) assumed an evolution of density from z=5.7 to z~7. They found that 



the IGM is not highly neutral at z~7 an d the neutra l hydr ogen fraction could evolve from 
x m 6 ' 6 ~ 0.24-0.36 to ar^/ ~ 0.32-0.64. iHibon et all (boilh has also performed a search for 
z~7 LAEs using the IMACS instrument at the Magellan telescope. This study resulted in a 
sample of 6 z~7 LAEs candidates, for which the spectroscopic follow-up data are in analysis. 



We present here a new NB imaging survey with the Suprime-Cam/Subaru telescope 
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- targeting z ~ 7 Lyman-a emitters with the custom-made NB973 filter (A cen ter=9755A, 
FWHM=200A). This paper first presents the data (Section 1) and the data reduction proce- 
dure (Section 2). We then describe the method of selection and contamination of low-redshift 
interlopers for high redshift LAEs in Section 3. We present the final sample of z ~ 7 LAEs 
and Lya luminosity function at this redshift in Section 4. 



Throughout this study, we adopt the following cosmological parameters : H 
= 0.3, Q A = 0.7 f lSpereel et all 120070 . All magnitudes are AB magnitudes. 



70km. s l .Mpc 



2. Observations and Data Reduction 
2.1. Observations 

The data were taken with the Suprime-Cam instrument (the Subaru Prime Focus Cam- 
era), installed at the 8m Subaru telescope at the National Astronomy Observatory of Japan. 
This instrument delivers a mosaic of ten 2048 x 4096 CCDs, which covers a 34' x 27' field 
of view with a pixel scale of 0.20". The new CCD, with a QE two times better than the 
precedent one at the red end (A = 8000 - 10400A), was installed. 

By taking exposures lasting 20 min, we ensured background limited performance. We did 
offset the telescope between each exposure. The offset sequence, for each field and night, 
was randomly chosen. 

We targeted two fields (hereafter called D33 and D41) covered by the Canada France Hawaii 
Telescope Legacy Survey (CFHT-LS) and the WIRCam Deep Survey (WIRDS : Pis Willott 
& KneibjE The total area of the survey, the two fields considered, is 2340 square arc minutes. 

We observed using the custom-made NB filter centered at 9755A with a wavelength 
range of AA=200A (iV\B973) during two nights in July 2009. During these observations, the 
conditions were partially good : we obtained 15 and 19 20-minute exposures during these 
nights with a seeing varying between 0.5" and 0.85" . 
The total exposure time for each epoch of data is given in Table [TJ 



1 Based on observations obtained with WIRCam, a joint project of CFHT,Tai wan , Korea, Canada, France, 
at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Research Council (NRC) 
of Canada, the Institute National des Sciences de l'Univers of the Centre National de la Recherche Scientifique 
of France, and the University of Hawaii. This work is based in part on data products produced at TERAPIX, 
the WIRDS (WIRcam Deep Survey) consortium, and the Canadian Astronomy Data Centre. This research 
was supported by a grant from the Agence Nationale de la Recherche ANR-07-BLAN-0228 
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2.2. NB Data Reduction 



Using the Data Reduction Softw are developed for the Subaru Suprime-Cam instrument 
( lOuchi et al.ll2004j ; lYagi et al.l 120021 ) . we performed the following steps: a bias subtraction, 
a flat fielding, a distortion and atmospheric dispersion correction, a sky subtraction, a bad 
regions masking such as satellite trails and AG probe, the alignment of the individual expo- 
sures, and the final co-adding step resulting on a stacked image. After eliminating the low 
signal-to-noise regions at the edges of the field of view, we obtain an effective area of 1118 
arcmin 2 for D33 field, and 1202 arcmin 2 for D41 field. 

We realize the astrometr ic calibration on the i ndividual images before the final stacking us- 
ing the UCAC2 catalog (IZacharias et all 120041 ) . 

We need then to adjust the WCS parameters to obtain a more precise alignment. For this 
purpose we use the IRAF task msctpeak to interactively align the catalog stars in our images 
and update their headers. We finally obtain an astrometry calibration for each individual 
images with a precision of rms~ ±0.1 arcsec in both directions. 

The photometric calibration of the u*, g', r', i', z' broad band CFHT-LS data (described in 
the next section) is based on the SPSS data for stars with 1 7< i' <21 and the Megacam- 
SDSS color transformation equations of iRegnault et al.l ( 120091 ) .The precision obtained in u*, 
g', r', i', z' is between 0.03 and 0.02 mag. As the NB973 filter is included in the ^'-band 
filter, we calibrated our NB973 stacked images, using the AUTO magnitude from the z' 
band SExtractor catalog. We performed the calibration on 1500 non-saturated stars within 
16< z' <20. Considering the photometric error on the broad-band calibration, we obtain a 
photometric calibration precise at 0.1 magnitude in NB973. 

To estimate the limiting magnitude for the different bands, we used the script limitmag.cl 
available in the SDFRED package. This task constructs a count distribution from random 
photometric apertures on the image. Then it fits a Gaussian profile into this distribution, 
obtains the sigma number of this profile and calculates the limiting magnitude of the image. 
We report in Table [1] the limiting magnitude for each band in each field. 



2.3. Broad Band Data 

Very deep optical imaging data of our observed fields are available through the CFHT- 
LS. For the purpose of this study, we made use of the T0006 release. These data products 
are available from the CADC archive and take form of image stacks in the u*, g', r', i', z' 
filters and of ancillary data such as weight maps, catalogs etc. The spectral curves of the 
filter u*, g', r', i', z' are similar to the ones of SDSS filters. We also have in hand the Near 



-5 - 



Infrared deep Imaging data from the WIRCam Deep Survey (WIRDS) (Pis: CWillott, J-P. 
Kneib) (Bielby et al. in preparation). 

These optical data have been calibrated photometrically using the SDSS photometry and 
the NIR data using 2MASS photometry (IMcCracken et al.l 120101 ) . Considering internal and 
external photometric error sources, the uncertainty on the optical and the NIR data pho- 
tometry is ~0.03 mag and ~0.02 mag, respectively. 

As the broad-band data and the narrow-band data do not have the same pixel scale, we re- 
sampled the broad-band data using the software SWARP to obtain optical and NIR images 
with 0.2 arcsec/pixel. The alignment in pixels is then verified using IRAF geomap/geotran 
tasks. 

Our complete set of data is therefore scaled at 0.2 arcsec/pixel and is covering an effective 
area of 1118 and 1202 arcminutes 2 , for D33 and D41 fields, respectively. 



A summary of the observational data for both fields, D33 and D41, used in this paper 
is provided in Table [TJ Figure [1] shows the transmission curves of the filters corresponding 
to the multi-band data used in this study. 



3. Sample 
3.1. Catalog generation 



We generate the catalogs using the software SExtractor of iBertin fc Arnoutsl ( 119961 ). 
We used the dual-image mode : the first image, the detection image, has been settled as 
the combined NB973 image, the second image, the measurement image, corresponds to the 
resampled images from the optical and NIR bands. 

We choose to detect objects in 7 pixels above a threshold of 1.2er. The aperture used for the 
photometry is 1". 



3.2. Selection 

Criterion#l : SNR(A^B973 combined )> 5a. 

We selected objects with a 5a detection on the combined NB973 images. 



Criterion#2 : SNRfV, g' , r', i', X 2 )< 3d. 

Due to the Gunn- Peterson trough, i.e. the nearly complete absorption of the flux short- 
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Fig. 1. — Transmission curves of the filters corresponding to the complete set of data used 
in this paper. All transmissions are normalized to 100% at maximum. 
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ward of Lymana as a result of the large neutral hydrogen column density in the Intergalactic 
Medium (IGM), we observe spectral discontinuity at redshifts greater than about 6. We are 
therefore searching for objects which are not detectable in optical (u*, g', r', i') bands. A 
possible method is therefore to select objects with less than a 3a detection in filters blue- 
ward of the expected Lya emission : u*, g', r', %' . 

The color break between the optical and NB973 filters is high and covers a wide spectral 
range. Moreover, for the CFHT-LS, the Terapix data center generated deep x 2 image com- 
bining the g', r' and i! images. We consider therefore this x 2 image as well as the other 
blue-ward optical filters and we select objects with less than a 3a detection in the x 2 image. 
In summary, we selected objects with less than a 3a detection in u*, g', r', z',and x 2 bands. 



Criterion#3 : a color criterion between z' and NB973 data. 

The NB973 filter used for this study is included in the broad z band filter. Although we are 
expecting an excess of flux in NB973, it is yet possible to observe part of the continuum in the 
z band filter. From Criterion^! and Criterion^, we can define the following col or criterion 
: z' - NB973 > O.Q5mag. We used the same NB973 filter than bta et~aP koO$ \ have used 
for the search of robust z=7 LAEs in the SDF field. Their search was successful as it resulted 
on the spectroscopic confirmation of the only z~6.96 spectroscopically confirmed LAE. Fol- 
lowing the simulation they performed and plot on their Figure 3, we also check the color 
z'— NB973 of our candidates with a second color criterion : z' — NB973 > 1.72mag. We show 
our candidates in the Figure [2] representing z' — NB973 vs NB973 for both fields, and the two 
color criteria. The points answering z' — NB973 > O.Q5mag but not z' — NB973 > 1.72mag 
are the points with no z' band detection. We used then the 3a limit of the z' band data for 
these objects. 



Criterion#4 : NB973 - J < 0, ^5973 - H < and NB973 - Ks < 0. 

To avoid some contamination by low-redshift objects, we choose to set NIR criteria. More 

explanations are given in the section 13.31 

Criterion#5 : High redshift objects are identified as compact point-like sources. After apply- 
ing carefully all the criteria presented above, we inspect each candidate visually to confirm 
their point-like source aspect. Objects with a different aspect are more likely artefacts. 7 ob- 
jects out of 9 for D33 and 7 objects out of 10 for D41, have been selected as serious candidates. 
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Table 1: Observational data. 







Integ 


ration 


Limiting 


Instrument 


Band 


time (hrs) 


magnitudes 3. 






D33 


D41 


D33 


D41 


MegaCam 


u* 


21.3 


21.4 


27.5 


27.4 


MegaCam 


g' 


22.1 


24.4 


27.5 


27.2 


MegaCam 




39.6 


40.7 


27.0 


26.8 


MegaCam 


a 


69.3 


65 


26.6 


26.5 


MegaCam 


z' 


48.7 


44.5 


25.5 


25.4 


Suprime-Cam 


NB973 


5 


6.3 


24.3 


24.7 


WIRCam 


J 


4.6 


4.8 


24.5 


24.4 


WIRCam 


H 


4.4 


4.1 


24.4 


24.1 


WIRCam 


K s 


4.8 


3.9 


24.2 


24.3 



5cr magnitude limits in apertures 2" in diameter for MegaCam, WIRCam and Suprime- 
Cam. 




22 22.5 23 23.5 24 24.5 22 23 24 25 



NB973 mag (AB) NB973 mag (AB) 

Fig. 2.— Color-Magnitude Diagram z' - NB973 vs 7V5973 for D33 field (left) and D41 
(right), showing the candidates (red triangles) obtained with the criteria z' — NB973 > 
0.65mag (green dashed line). We also compare them to the criterion z' — NB973 > 1.72mag 
(blue dashed line). The black dots represent all the objects present in each field with 
SNR(NB973 combin ed)> 5a. 



-9 - 



3.3. Contaminants 



3.3.1. Transient objects 



Since the narrow band and broad band data were taken at different times, it is possible 
for transient objects to appear in one filter and not another. Transient objects that are 
brighter than the narrow band detection limit will be considered as candidates if they are 
not visible in the broad band images. 

Kulkarni fc Raul ( 120061 ) used a transient rate fo r SNe, including Type la, b,c and Type II, of 
5 x 10 4 Gpc 3 yr _1 . From ICappellaro et al.l ( 19991 ). we know that the transient rate of Type la 
SNe is approximately a factor 3.4 lower than the transient rate including Type la, b, c and 
Type II SNe. Therefore, by applying a transient rate of 5 x 10 4 /3 to each of our fields, we 
can estimate that 0.8 Type la SNe could be detected in the D33 field and 0.9 Type la SNe 
in the D41 field. 



We chose to also estimate the number of SNe that we could have detected in our sur- 



vey using t he m ethod of lOta et al.l (120081 ) . We assume the same mean object colour than 
Ota et al.l (120081 ). i' — NB973 = .33, and u s e the results of variability. For their large 
magnitude variation (1.1-1.6 mag), lOta et al.l ( 120081 ) find one object at each epoch giving 
P ~ 10~5. Down to i' ~ 26.5, we found approximately 7500 objects in D33 and 9000 objects 
in D41. From this method, we obtain therefore an estimated number of variable objects of 
~ 0.75 object for D33 and ~ 0.9 object for D41 . 



In summary, both methods seem to agree for a possible contamination of order one 
transient object per field. 



3.3.2. L-T dwarfs 



Considering Figure 9 of lHawley et al.l (120021 ) presenting several color diagrams for L-, 
T- dwarfs detected in SDSS, we know that M-, L-, T- dwarfs have z' — J > 1 and L-, T- 
dwarfs have z' — J > 2. Therefore, although our J-band is shallower than our NB973 limiting 
magnitude, the criterion NB973 — J <0 still secure a non- contamination of our LAEs sample 
by L-, T-dwarfs. We therefore avoid selecting L- and T- dwarfs. 
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3.3.3. Foreground emitters 

We estimated the lower value on equivalent width a line emitter, wit h a flat continuum in 
'v, wo uld require to be selected with our criteria using the formula from iRhoads fc Malhotra 



( 1200 If ): 



EW„ 



NB 



BB 



AX 



NB 



5a 



NB 



3oBB 



AA 



NB 



(1) 



with and Jbb the flux in NB973 and g' band respectively, AAjvb the width of the NB973 
filter, o~nb and o~bb, the flux uncertainties in NB973 and g' band respectively. We obtained 
therefore an EW®f£ ~ 840A and EW^^ ~ 714. 2A. Foreground line emitters would then 
require an equivalent width EW^j^ > 840 A and EW®^ > 714.2 A to contaminate our Lya 
selection sample. 

1- Ha at z~0.47 



Salzer et al.l (120051 ) presents the typical equivalent width of low redshift galaxies. They 
found that Ha emitters at z~0.4 have EW < 200A. From Equation fTJ the lower limit on 
the Ha emitters equivalent width in our survey is 840A in D33 field and 714. 2A in D41 field. 
Ha emitters at z~0.47 can therefore not contaminate our candidates selection. 



2- [Om] at z~0.95 



Xia et al.l (120101 ) studied the line emitters in the PEARS survey and observed that [O ill] 



emitters at z~l have generally an equivalent width lower than 500A. Moreover, looking at 
the typical Spectral Energy Distributions (SEDs) of [O ill] emitters at z~0.9, we remark that 
these foreground emitters are brighter in Ks band than in J band and with typical J band 
magnitude ~ 23 (AB). We will have therefore observed these objects in our very deep NIR 
data. Considering the lower limit on equivalent width that line emitters will have to con- 
taminate our candidate sample and the NIR magnitudes of typical z~0.9-l [Om] emitters, 
the contamination by th ese foreground emitters is therefore ruled out. 
Cardamone et al.l (120091 ) found low redshift [O ill] emitters with very high observed EWs, 
reaching values until 1500A. However, these emitters have very bright magnitudes in optical 
bands : 18 < g',r',i',z' < 21 and blue spectra. They would be therefore easily detected by 



the CFHT-LS data. These emitters are identified as Pea galaxies 
galaxies. 



luminous blue compact 



- 11 - 



3- [Oil] at z~1.6 



Hayashi et al.l ( 1201 ll ) studied [O il] emitters at z~1.46. Figure 8 of lHayashi et al.l (j201lh 
shows these emitters in a colour magnitude diagram z' — K s vs K s . z~1.46 [O ll] emitters 
have K s e[20; 23], and z' — K s > 0. As our narrow-band filter is included in the z' band filter, 
as seen in Figure HJ we ca n assume that z' — K s > 0. This colour does not agree with our 
Criterion 7^4. Figure 10 of lHayashi et al.l (120111 ) shows the colour-colour diagram z' — K s vs 
J — K s . Most of their [O ll] emitters have not only a z' — K s > but also a J — K s > 0. In 
the case of our sample, a deeper J band image would help us eliminate the [O ll] emitters 
as a possible contaminant to verify that our candidates have J — K s < 0, as deduced from 
Criterion$4. 

However, our K s band data do not cover entirely both fields. Although the contamination by 
z~1.6 [On] emitters is very unlikely due to the K s magnitude range, we cannot completely 
rule it out. 



3.3.4- R e d Continuum Galaxies 



Balmer break galaxies 

Balmer break galaxies at z~1.4 need t o be also considered a s a possible contamination 
source. However, from Figure 4 from iDelorme et al.l (j2010[ ). we o btain an estimation 
of the colour z' — J for z~2 Balmer break galaxies : z' — J > 0. iDaddi et al.l (120041 ) 
shows that z> 1.4 Balmer Break galaxies have z' — Ks > 1.5. Moreover, in the case of 
continuum objects, we can assume that z' — Ks ~ NB973 — Ks. During the selection 
of our candidate sample, we applied Criterion$4 : z' — J < 0. If the candidates were 
Balmer break galaxies, we should have detected them in J with a brighter magnitude 
in J than in z' . 



Extremely Red Objects (EROs) 



We add a criterion to avoid selecting extremely red objects ( ICimatti et al.l 120021 ) . 
We constrain this contamination by applying NB973 — J < 0, and also verifying 
NB973 - H < and ^£973 - Ks < 0. 



These red continuum galaxies are therefore generally brighter in Ks band than in NB973. 
These objects cannot contaminate our candidate samples. However, as our NIR data do not 
cover our entire fields, we cannot completely rule out this contamination. 
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3.3.5. False Detections 

In order to estimate the number of false detections that could pass our selection criteria, 
we create an inverse NB973 combined image by multiplying this NB973 image by -1. We 
then applied the first criteria used for the candidate selection : SNR(iV.B973 com fcj nec ;)> 5a. 
None detection meets this condition. We found out therefore that our candidate samples are 
not contaminated by false detections. 



3.4. Summary 

After this analysis, we can conclude that it is unlikely that our z~7 LAE candidates is 
contaminated by low-redshift interlopers. 

A more detailed study of the foreground emitters (Ha at ~0.47, [Om] at ~ 0.95 and [On] 
at ~ 1.6) will be presented in a forthcoming paper. 



3.5. Final sample 

We obtain a final sample in each observed field. The D33 final sample contains 7 
z~7 LAEs candidates over the range of ^5973=24.3-24.8 and SNR (JV.B973)=5.2-7.7, and 
the D41 finale sample, 7 z~7 LAEs candidates in a range of iV\B973=24. 6-25.0 and SNR 
(iV\B973)=5. 6-7.9, as described in Table [2] and Table [3j Using the equation given below, 
we estimate the lower limits of the rest-frame equivalent widths (EW) derived from the 
photometric data and report them in Table |2] and Table |3j It is interesting to observe that 
the candidates from the D41 field, deeper than the D33 field (see Table [1]) are automatically 
fainter than the D33 ones. The observational conditions for both fields show a systematic 
difference of at least 0.1 arcsec in seeing. 



EW, 



rest 



fNB973A\ z > — f z 'A\ NB973 



fz> - f- 



NB973 



X 



l + Z 



(2) 



where fNB973 is the observed flux in the narrow-band combined image, f z > is the observed 
flux in the z' broad-band image, AAatb973 and A\ z > are the width of the NB973 filter (200A) 
and the z' band filter (928A) respectively. 



Following lOta et al.l (120081 ) , we assume that 77% of the NB flux comes from the Lya 
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line : fiy a ~ 0.77/^5973. We use the detection limit in z' band, to derive in turn lower EW 
limits. 

In the D33 field, LAE^l has a NB973 magnitude corresponding to 94% completeness 
and the LAE#7 ^5973 magnitude to 89% completeness. In the D41 field, LAE#1 has 
a NB973 magnitude corresponding to 96% completeness and the LAE^7 NB973 magni- 
tude to ~89% completeness. To estimate the completeness, we added 200 artificial star-like 
objects per bin of 0.1 magnitude in blank regions of the stacked images. We then run SEx- 
tractor on the image with the same parameters as previously used for object detection. This 
procedure has been repeated 20 times. The average count on 20 times of the number of artifi- 
cial stars retrieved in each magnitude bin provides a direct measure of the completeness limit. 



Discussion 



4.1. Variance 



Two sources of variance can affect a high redshift study : the Poisson variance and the 
cosmic variance, in other words the fluctuations in the large scale distribution of the galaxies. 
In order to e stimate the value o f this cosmic variance we used the on-line calculatoi@ from 
the model of iTrenti fc StiaveUil (120081 ). We obtain then a value of 37% and 36% of cosmic 
variance in D33 and D41 fields respectively. 



2 http://casa. colorado.edu/ trenti/Cosmic Variance. html 



Table 2: Table of the z~ 6.96 LAE candidates for the D33 field. 



Name 


^973 


Error 


SNR (^5973) 


z' 


Error 


SNR (z') 


EW 3 - (A) 


LAE#1 


24.3 


0.17 


7.7 


> 25.46 






> 24 


LAE#2 


24.6 


0.2 


6.1 


> 25.46 






> 24.3 


LAE#3 


24.7 


0.21 


5.7 


> 25.46 






> 24.3 


LAE#4 


24.7 


0.21 


5.6 


> 25.46 






> 24.3 


LAE#5 


24.7 


0.22 


5.5 


:> 25.46 






> 24.3 


LAE#6 


24.8 


0.22 


5.3 


26.5 


0.31 


3.7 


1 


LAE#7 


24.8 


0.23 


5.2 


26.8 


0.38 


2.8 


7.6 



In the rest-frame 
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LAE#1 LAE#2 LAE#3 LAE#4 LAE#5 LAE#6 LAE#7 




Fig. 3. — Thumbnail images of the D33 candidates listed in Table EJ Each window is 
15" x 15". Objects names and passbands are located above and to the left of the thumbnails, 
respectively. 




Fig. 4. — Thumbnail images of the D41 candidates listed in Table EJ Each window is 
15" x 15". Objects names and passbands are located above and to the left of the thumbnails, 
respectively. 



- 16 - 



Considering the limited number of objects in our sample and the large comoving volume 
of our survey, our results are as limited by Poisson noise - ~38% for the 7 objects of D33 
field and for the 7 objects of D41 field - than by clustering. Therefore, in the case where 
the full sample in each observed field (7 candidates in each field) is taken into account (see 
Par agr aph 14 . 2 . 1 p . the error bars represent the Poisson variance and the cosmic variance. In 
the two other cases (see Paragraphs 14.2.21 and I4.2.3[) . our results are more limited by the 
Poisson noise than by clustering. The error bars are therefore representing the Poisson noise 
only. 

For D33 field, we obtain a total fractional error on number counts of 0.65. As we obtain 
7 objects in D33, we should expect between 4.5 and 9.4 objects in D41, if the field-to-field 
variation is within la. As we found 7 objects D41 field, we can therefore confirm that the 
field-to-field variation is within la. This result was expected as none of our field are known 
to be located in overdense or underdense large-scale structure at this redshift. 



4.2. Luminosity Function 



Following lOta et al.l (120081 ). we assume that on average 77% of the narrow-band flux 
comes from the Lya line. We therefore apply the same correction factor during the conver- 
sion from NB973 to Lya fluxes. 



We fit to the Lya luminosity function of these z~7 LAE samples, a Schechter function, 
$(L), given by 



Table 3: Table of the z~ 6.96 LAE candidates for the D41 field. 



Name 


^5973 


Error 


SNR (^5973) 


z' Error 


SNR(z') 


EW A (A) 


LAE#1 


24.6 


0.17 


7.9 


> 25.4 




> 24.3 


LAE#2 


24.7 


0.18 


7.1 


> 25.4 




> 24.4 


LAE#3 


24.8 


0.19 


6.6 


> 25.4 




> 24.5 


LAE#4 


24.9 


0.2 


6.0 


> 25.4 




> 24.5 


LAE#5 


24.9 


0.3 


5.9 


> 25.4 




> 24.5 


LAE#6 


24.9 


0.21 


5.8 


26.7 0.38 


2.8 


4.4 


LAE#7 


25.0 


0.22 


5.6 


> 25.4 




> 18.4 



In the rest-frame 
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in order to compare w ith previous high redshift works ( IHibon et al.ll2010l ; lOuchi et al.ll2010 



2008; Ota et al 



20081 ). The error bars shown in Figure |5] represent the total errors (cumu- 



lative Poisson errors and cosmic variance) as explained in the previous section 14.11 In the 
case of Figure Eland the results are more limited by the Poisson noise than by the cosmic 
variance. The error bars represent then only the cumulative Poisson errors. Considering the 
low number of candidates in o ur sample, we choos e to fit two out o f three of the Schechter 
function parameters. Following IHibon et al.l (120101 ) and lOuchi et al.l ( 120081 ) . we set the faint 
end slope of the luminosity, a, to a = —1.5, and determine $* and L* by \ 2 minimization. 
We corrected the data from the completeness by number weighting. 



We decided to study the following different cases : 



• The full sample of each field is real. 

• The samples are contaminated by 50% . We then consider only the 4 brightest objects 
of each field as real (referred later as bright samples). This is justified by the fact that 
the brightest candidates are the most robust ones and we are focusing on building the 
bright end of the z~7 Lya luminosity function. 

• We chose to derive a common Lya LF for the D33 and D41 LAEs brightest candi- 
dates. We have therefore a sample of 8 bright objects in a total area of 2320arcmin 2 , 
corresponding to the sum of the D33 and D41 effective areas. 



The best-fit Schechter LF parameters for each case are summarized in Table HI 



4-2.1. All the candidates are real. 



From the z~7 Lya LF derived from the full sample of the D41 field, seen in Figure [5] 
as the blue line and the candidates corresponding represented as filled s quares, we ob s erve a 
possible but weak evolution in luminosity from the z~6.5 Lya LF from lOuchi et al.l (120 10f ) 
(dot-long dashed line in Figure |5J), but no possible evolution from the z~6.5 Lyct LF from 
Kashikawa et al.l (120111 ) (dot-short dashed line in Figure E]). The z~7 Lyct LF from the 
complete sample of the D33 field, shown as the red line and the filled triangles in Figure [5] 
shows a weaker evolution in luminosity from the z~6.5 Lya LF. However, looking at the 
error bars from D33 and D41 field candidates, both weak evolutions in luminosity, become 
very discussable. From these possible z~7 Lya LFs, if all the candidates are real for one 
of the observed field, we cannot therefore conclude about a possible evolution in luminosity 
from z~6.5 to z~7. 
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Fig. 5. — Best-fit Schechter function for the cumulative z~7 Lya luminosity function de- 
rived from the D33 full sample in red and from the D41 full sample in blue. Our D33 
candi dates are repres ented as triangles and the D41 ones as fille d squares, the z=6.96 LAE 
from live et all (120061 ) is the empty square and lOta et all (120081) candida tes are pentagons. 
Also represented here is the z =6.5 Lyq LF from iKashikawa et all (120111 ) (dot-short dashed 
line), th e z=6.5 Lya LF f r omlOuchi et all ( 120101 ) (dot-long dashed line) and the z=5.7 Lyct 
LF from IKashikawa et all ( 120111 ) (dotted line). The cumulative Poisson errors and the cosmic 
variance are taken in account in the vertical error bars. 
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Table 4: Best-fit Schechter LF parameters for a = —1.5 



Redshift 


log(L*(ergs 1 )) 


log($*(Mpc" 3 )) 


6.96 (1) 


42.50i°^ 


-2.40+°; 2 5 


6.96® 


42.53l°i 


-2-151°;} 


6.96® 


42.731^ 


o i 7+0.15 
O.Li _ 2 


6.96® 


42.56+°;! 


-2.731°-} 


6.96® 


43.71±^;i 


-3.33«;} 


6.96® 
6.5® 


42 8+°- 12 
42.64+°-} 


_o 44+0.20 

_o Q7+0.13 
°- u ' -0.13 


5.7® 


42.8±g:ig 


q i i +0.29 
°- ii -0.31 



References. (1, 2) derived from our D33 and D41 full samples resp.; (3, 4) derived from the 
brightest can didates in each fie ld; (5 ) derived from a co mmo n sample between D 33 and 
D41 fields; (6 jota et al.l J2OO8I); (7) lOuchi et al.l ( hoioh : (8) lOuchi et al.l ( 120081 ) 



4-2.2. Only the brightest candidates are real. 



We derived Lya luminosity functions for the two different z~7 bright LAE candidate 
samples. The best-fit Schechter functions are shown in Figure [6j In both cases, we do 
not observe a sign ific ant evolution either in density nor in luminosity from z~6.5, from 



Ouchi et all (120101 1 or lKashikawa et al.l (120111 ). to z~7. 



4-2.3. One common LF for both fields bright sample. 



By fitting a unique Schechter LF for both samples (as seen in Figure [7]), we assume 
that the brightest objects from both samples are real z~7 LAEs. Looking only at the filled 
triangles in Figure [7J representing the four brightest candidates from D41 and the o nes from 
D33, we remarked that t hey are fitting the z~6.5 Lya LF from lOuchi et al.l (120101 ) but not 
the z~6.5 Lya LF from iKashikawa et al.l ( 120111 ) . Depending on the z~6.5 Lya LF we 
considered, the best-fit Schechter derived LF, seen as the red line on Figure El could agree 
with the observed evolution between z~5.7 and z~6.5 but would not be in favor of a strong 
evolution between z~6.5 and z~7. 
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L(Lya) (erg s" 1 ) 



Fig. 6. — Best-fit Schechter function for the cumulative z~7 Lya luminosity function derived 
from the D33 bright sample in red and from the D4 1 bright sample in blue. Our D33 
cand idates are repres ented as triangles and the D41 ones as filled squares, the z=6.96 LAE 
fr om live et al.l (120061) is the empty square. Also represented here is the z=6.5 Lyq LF from 



Kashikawa et al. 



(120 111 ) (dot-short-dashed line), the z=6.5 Lya LF fromlOuchi et al.l (120 lOh 
(dot-long-dashed line) and the z=5.7 Lya LF from iKashikawa et al.l (j201l[ ) (dotted line). 
The vertical error bars represent the cumulative Poisson errors, as the bright samples are 
more limited by the Poisson noise than by clustering. 
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L(Lya) (erg s" 1 ) 



Fig. 7. — Best-fit Schechter function for the cumulative z~7 Lya luminosity function derived 
from a common sample composed by the D33 bright sample in blue and from the D41 bright 
samp le in red. Our D33 and D41 candidates are represented as triangles, the z=6.96 LAE 
from live et al.1 d2006h is the empty square. Also represented here is the z=6.5 Lyq LF from 



Kashikawa et al. 



(1201 If ) (dot-short-dashed line), the z=6.5 Lya LF fromlOuchi et all (120101 1 
(dot-long-dashed line) and the z=5.7 Lya LF from iKashikawa et al.l (j201lh (dotted line). 
The vertical error bars represent the cumulative Poisson errors, as the common sample is 
more limited by the Poisson noise than by clustering. 
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4-2.4- Summary 



We obtain two different best-fit Schechter functions for the z~7 cumulative luminosity 
function from our two different sets of photometric candidates. Both candidate samples help 
for building the bright end of the z~7 LF. Although we do not see an evolution from z~6.5 
and z~7 from the D33 candidates sample, the LF derived from the D41 candidates sample 
do not show a significant evolution from z~6.5 to z~7, as seen in Figure |6j 
By producing z~7 LF including not only a one-field candidate sample but also the spectro- 
scopically confirmed IOK-1, the best-fit Schechter parameters do slightly change. However 
our previous conclusion is still valid. 

Assuming that the 8 brightest of our candida tes are real, we int erestingly obtain a z~7 LF 
in agreement with the z=6.5 LF produced by lOuchi et al.l (120101 ) . 



Ota et al.l ( 120081 ) presents the first search for z~7 Lya emitters which has lead to the 



first spectroscopically confirmed z~7 LAE. Their imaging survey covers an area of 876 square 
arcmin with the filter NB973 (AA = 200A, A c = 9755A) and reaches a 50% completeness 
of ^5973=26.2 (AB, 5a). IOK-1 has a flux of 2 x 10~ 17 erg s" 1 cm" 2 . We therefore observe 
a wider but shallower area. This strategy is justified by our willing to constrain only the 
bright-end of the z~7 Lya LF. 



The main prediction of iLe Delliou et al.l (120061 ) is a likely moderate decline of the bright 
end of the LF of LAEs from z=6.5 to z~7 arising from the evolution of the mass distribution 
of dark matter halos. If the LF may undergo limited evolution between z=6.5 and z=7, the 
effects induced by the incomplete reionization of the IGM may play an important role in 
the evolution of the observed LF. The faint galaxies may be more easily obscured by neutral 
regions and coul d enhance the bright end of the LF. 

Ota et al.l (120081 ) . by assuming that the LAEs and LBGs show a similar evolution history, 
can infer that the L ya LF derived from LAEs evolves as the rest-frame UVLFs obtained 
from LBGs. From lYoshida et al.l (120061 ). they estimated therefore a possible z~7 Lya 
LF with a pure lumino sity evolution of L* z=7 = 0.58L* =5 7 , with L* z=5 7 = 1.08.10 43 erg s _1 
( IShimasaku et al.l 120061 ) . This inferred z~7 Lya LF does not agree with our photometric 
candidates sample. It confirms their idea that, if LAEs are strongly related to LBGs, the 
neutral hydrogen fraction in the IGM is possibly higher at z~7 than at z~5.7, and this 
difference can cause the attenuation of Lya lines of high redshift LAEs. 



If none of our object is a real z~7 LAE, we can then put an upper limit on the z~7 
Lya LF, which helps for constraining better these models. 
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5. Conclusions 

We observed 0.64 square degree of the WIRDS/CFHT-LS fields with the NB973 filter, 
corresponding to targeting the Lya line at z~7. After applying our selection criteria and 
verifying that our selection was not contaminated by low-redshift emitters, we obtained a 
sample of seven z~7 LAEs in each field. From these photometric samples, we have been 
able to infer possible z~7 Lya Luminosity Functions. We did not find a significant evolution 
either in luminosity nor in density from z=6.5 to z~7. 

It is now crucial to obtain spectroscopic follow-up observations to reveal the real nature of 
these objects and establish a firm conclusion on the z~7 Lya Luminosity Function. The 
exact evolution of the LF beyond redshift 6.5 remains therefore a matter to debate. More- 
over, a single galaxy at z~7 is obviously not sufficient to firmly constrain the bright end of 
the Lya LF. We need to increase the number of independent fields as well as the number 
of LAEs at z—7. Once the bright end of z~7 Lya LF is determined and possible evolution 
from z=6.5 is derived, it will become easier to assess whether the 1.06/zm and the 1.19/xm 
NB filters could reveal z~8-9 LAEs. 
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